The interlayer cohesive energy of graphite from 
thermal desorption of polyaromatic hydrocarbons 
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We have studied the interaction of polyaromatic hydrocarbons (PAHs) with the basal plane 
of graphite using thermal desorption spectroscopy. Desorption kinetics of benzene, naphthalene, 
coronene and ovalene at sub-monolayer coverages yield activation energies of 0.50 eV, 0.85 eV, 1.40 eV 
and 2.1 eV, respectively. Benzene and naphthalene follow simple first order desorption kinetics while 
coronene and ovalene exhibit fractional order kinetics owing to the stability of 2-D adsorbate islands 
up to the desorption temperature. Pre-exponential frequency factors are found to be in the range 
10 14 -10 21 s _1 as obtained from both Falconer-Madix (isothermal desorption) analysis and Antoine's 
fit to vapour pressure data. The resulting binding energy per carbon atom of the PAH is 52±5 meV 
and can be identified with the interlayer cohesive energy of graphite. The resulting cleavage energy 
of graphite is 61 ± 5 meV/atom which is considerably larger than previously reported experimental 
values. 

PACS numbers: 81.05.Uw,61.50.Lt,68.43.Mn,81.07.De 
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I. INTRODUCTION 

The cohesive energy of a solid is determined by the 
interactions among its constituent atoms. In the case 
of graphite, a layered material with strongly anisotropic 
bonding it appears natural to distinguish the total from 
the interlayer cohesive energy. The first is dominated by 
strong localized covalent bonding through sp 2 carbon or- 
bitals while the latter is dominated by weak non-local van 
der Waals (vdW) interactions between graphene sheets. 
A description of the cohesive energy of graphite thus nec- 
essarily involves interactions of fundamentally different 
character and poses a true challenge even to the most ad- 
vanced calculational techniques. In particular dispersion 
forces that give rise to the long-range attraction between 
graphene layers have been notoriously difficult to predict. 
Not too surprisingly, one finds that values calculated 
from semi-empirical or ab inito methods for the interlayer 
cohesive or exfoliation energy of graphite range from as 
little as 8mcV/atom up to 170 mcV/atoimVi 3 : 4 : 5 : 6 Ex- 
perimental determinations of the interlayer cohesive en- 
ergy of graphite have been comparatively rare and are 
restricted to a heat of wetting experiment by Girifalco 
which yields an exfoliation energy of 43 meV/atomZ£ and 
a measurement by Benedict et al. based on radial de- 
formations of multi-wall carbon nanotubes which yields 
35meV/atomiS At this point it seems that not only the 
agreement between theory and experiment leaves room 
for improvement but that the experimental evidence for 
such comparisons should also be put on a firmer basis. 

Here, we aim at a better experimental characteriza- 
tion of the weak interlayer interactions in graphite using 
thermal desorption of thin films of polyaromatic hydro- 
carbons from the surface of a highly oriented pyrolytic 
graphite (HOPG) sample. Polyaromatic hydrocarbons 
(PAHs) are planar aromatic molecules formed by two 



or more fused aromatic rings and valencies of periph- 
eral atoms are satisfied through covalent bonded hydro- 
gen atoms. The interaction of PAHs with graphite is 
thus considered as suitable model system for the inter- 
action between graphene layers. This analogy is also 
supported by their striking structural resemblance, i.e. 
the same hybridization of atoms and practically identical 
bond lengths as well as the formation of adsorbate layers 
commensurate with the substratoi2iAiii£ which suggests 
that the character of the interaction between PAHs and 
graphite - in particular for larger PAHs - should be the 
same as that of the interaction of a graphene sheet with 
a graphite substrate. Moreover, the electronic structure 
and density of states of larger PAHs converges rapidly to 
that of graphene ~ 

A better understanding of the long range vdW forces 
in graphitic systems is also of interest when interac- 
tions between carbon nanotubes (CNTs) are studied, 
for example. The latter are frequently found to be ag- 
glomerated in so called carbon nanotube ropes, quasi- 
crystalline arrangements of close packed CNTs which are 
difficult so separate due to considerable long range vdW 
interactions^ Only recently have soaps been successfully 
used to separate and exfoliate such ropes^ where again, 
vdW forces between the surfactant and the tubes play a 
crucial role as for the wider field of soft matter physics. 



II. EXPERIMENTAL 

Thermal desorption (TD) experiments were performed 
under ultra-high-vacuum (UHV) conditions where the 
base pressure of below 1 x 10 -10 mbar was maintained 
using a combination of membrane, turbo-drag and turbo- 
molecular pumps. The HOPG sample from Advanced 
Ceramics (Grade ZYB) was mounted on a Ta disk using 
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conducting silver epoxy and was freshly cleaved prior to 
transfer into the vacuum chamber. The sample surface 
was cleaned prior to dosing by repeated annealing cycles 
to 1200 K. A type-K thermocouple was spot-welded to 
the Ta disk to measure and allow control of the sam- 
ple temperature. The thermocouple was calibrated using 
desorption of Xc multi-layers in combination with their 
well known heat of sublimation. The sample holder was 
attached to a He continuous flow cryostat that enabled 
sample cooling down to 30 K. 

Benzene or naphthalene (99.89% and 99.99%, Aldrich) 
vapour was admitted from a gas reservoir to the sample 
surface through a retractable pin-hole doser. Unwanted 
atmospheric contaminants were removed from the sol- 
vents by freeze-pump cycles prior to adsorption experi- 
ments. Exposure of HOPG to coronene (99%, Aldrich) 
and ovalene (99.5%, Dr. Ehrcnstorfer GmbH) was by 
sublimation of powder material from a PID temperature 
controlled Knudsen cell. A typical coverage series of TD 
traces was obtained by dosing the graphite surface with 
successively increasing adsorbatc quantities up to a to- 
tal coverage of approximately 5 ML. TD spectra were 
recorded for a constant heating rate between 0.75 K/s 
and 2 K/s. Desorption of species with a mass to charge 
ratio of up to 200 a.m.u/e, such as benzene and naph- 
thalene, could be monitored using a quadrupole mass 
spectrometer (Spectra Satellite). Due to high molecu- 
lar mass, coronene and ovalene desorption signals had 
to be obtained from the total pressure inside the UHV 
chamber as monitored by a Bayard-Alpert ionization 
gauge. Background correction of the latter TD traces 
was carefully cross-checked with simultaneously recorded 
TD-spectra from species with masses between 4 to 200 
a.m.u. Further details of the experimental procedure are 
available elsewhere^ 



III. RESULTS AND DISCUSSION 

A. Thermal desorption of PAHs 

We begin our discussion with benzene TD spectra 
shown in fig.^a,), which were recorded after exposure of 
the HOPG surface with a dose of up to 5 Langmuirs 
(1 L = 1 x 10~ 6 Torr s). The TD spectra are charac- 
terized by a high temperature feature attributed to des- 
orption from the first monolayer and a low temperature 
feature due to the desorption from multi-layers, which 
are each centered around 152 K and 142 K, respectively. 
The sub- monolayer phase diagram of benzene on graphite 
shows a complex pattern of phase transitions were in 
the high coverage- and high temperature- phases benzene 
molecules tend to be oriented with their molecular plane 
perpendicular to the graphitic surface to satisfy steric 
and entropic constraints *1 However, at coverages signif- 
icantly below one monolayer and at 152 K, i.e. the tem- 
perature of desorption in our experiments, the molecules 
are expected to be adsorbed with their aromatic rings 
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FIG. 1: Series of thermal desorption spectra of a) benzene, b) 
naphthalene, c) coronene and d) ovalene from HOPG. The 
rate of desorption was 0.7 Ks _1 , 1.0 Ks _1 , 2.0 Kb" 1 and 
2.0 Ks _1 respectively. Coverages range from the submono- 
layer regime to about 5 monolayers. 
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parallel to the surfaced The shape of TD spectra and 
the coverage independent desorption peak maximum at 
low coverages are clearly indicative of first order kinetics. 
With coverages approaching the monolayer regime, how- 
ever, the desorption traces broaden towards the low tem- 
perature side and become more complex - as expected 
due to the changes of molecular orientation within the 
adsorbatc layer at higher coverages. The TD spectrum 
from a complete monolayer is here associated with the 
last trace before the additional desorption feature around 
140 K develops with increasing coverage. The latter is 
typical of desorption from multi-layers with a common 
leading edge corresponding to zero order kinetics. 

As in the case of benzene, the coverage series for naph- 
thalene desorption exhibits two clearly distinguishable 
desorption features corresponding to mono- and multi- 
layer desorption respectively (see Fig.^i). The TD spec- 
tra of coronene and ovalene in Fig^:) and d) on the 
other hand do not exhibit clearly distinguishable mono- 
and multilayer desorption features. At low exposure, des- 
orption traces from both substances exhibit a behaviour 
indicative of fractional order kinetics. Here, the satura- 
tion of the first monolayer is estimated from the onset 
of multilayer desorption as apparent from the develop- 
ment of a common leading edge in the TD spectra at 
higher coverages. This leaves us with considerable un- 
certainty of the actual monolayer coverage which - as a 
worst case scenario - is assumed to be at most a factor 
of two higher or smaller. For ovalene, further evidence 
for the formation of the first monolayer is obtained from 
the appearance of a shoulder in the TD spectra before 
the common leading edge develops. The monolayer des- 
orption maxima for coronene and ovalene are found at 
390 K and 490 K, respectively. Adsorption of coronene 
monolayers on graphite reportedly also leads to the for- 
mation of 2-D islands^ owing to lateral interaction be- 
tween adsorbed molecules, which in the case of coronene 
and ovalene may make islands stable up to typical desorp- 
tion temperatures. Under such circumstances, diffusion 
from the island edges may become rate limiting and give 
rise to fractional order desorption kinetics as observed for 
coronene and ovalene (see analysis in the next section). 

B. Determination of frequency factors 
and binding energies 

The rate of desorption of some adsorbate from a solid 
surface is commonly described by an Arrhenius equation: 



— = -v6 n exp(-E a /k B T) (1) 
at 

where, v is the pre-exponential frequency factor, 9 is the 
surface coverage and E a is the activation energy for des- 
orption. 

In the following we will determine the activation en- 
ergies for desorption of each adsorbate using the series 



of TD traces presented above. Based on the above Ar- 
rhenius equation, a variety of techniques, each one with 
unique merits, can be used to obtain activation energies. 
Here, we will focus on an analysis using Redhead's peak 
maximum method^ and an isothermal analysis intro- 
duced by Falconer and Madi»i£. 

1. Redhead analysis 

The Redhead equation^ relates E a , the temperature 
at the desorption peak maximum T max and the heating 
rate, j3 as follows: 

E a = k B T max (ln(v T max //3) - 3.64) (2) 

The use of the latter expression is most commonly ap- 
plied to systems with first order kinetics but can be ex- 
tended to fractional- or zero order kinetics if the desorp- 
tion trace used for determination of T max corresponds to 
evaporation from a saturated monolayer. 

A crucial factor for the analysis of TD spectra using 
Redhead's peak maximum method is the availability of 
reliable pre-exponentials. Commonly, the latter are as- 
sumed to be on the order of 10 13 -10 15 s _1 and small un- 
certainties of less than an order of magnitude will not 
give rise to any serious error of the activation energy. 
However, the assumption of constant pre-exponentials is 
sometimes problematic j2i in particular if large adsorbates 
with many internal degrees of freedom are studied^i For 
thermal desorption of alkane chains from graphite for ex- 
ample, pre-exponential factors have been calculated us- 
ing transition state theory (TST) and were found to in- 
crease by over four orders of magnitude with chain length 
from 10 12 s _1 for methane up to 10 16 s~ 1 for Ci2H 2 6>2i 
Qualitatively this can be attributed to constrains on var- 
ious vibrational degrees of freedom of the molecule in 
the adsorbed state. If pre-exponentials are calculated 
from TST one generally computes the ratio of partition 
functions in the adsorbed and the transition state. If 
molecular degrees of freedom are somewhat constrained 
or 'frozen in' in the adsorbed state this will usually re- 
duce the corresponding partition function and thus tends 
to give higher pre-exponential factors. When studying 
and analyzing the thermal desorption of small and larger 
polyaromatic compounds we will thus avoid using esti- 
mated pre-exponentials and instead use experimentally 
determined values. 

We here assume that pre-exponential frequency fac- 
tors do not depend strongly on the film thickness and 
that multilayer values can thus also be used for analy- 
sis of monolayer desorption traces. Pre-exponentials for 
multilayer desorption can be obtained using the tempera- 
ture dependence of the adsorbate vapour pressure. This 
is done by assuming detailed balance between the rate 
of adsorption and desorption between a multilayer film 
in equilibrium with its gas-phase vapour. A similar ap- 
proach has previously been used by Schlichting and co- 
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workers The pre-cxponcntial frequency factor can then 
be expressed as: 

v = — = po (3) 

(T\/27r?TifcsT 

where s is the sticking coefficient, a the number of adsor- 
bates per unit area, m their mass and po is the vapour 
pressure at infinite temperature. The latter is obtained 
from the temperature-dependence of the vapour pressure 
using: 

p(T)=p exp(-^j (4) 

Here, AH S is the heat of sublimation at the tem- 
perature of desorption T max . However, none of the 
vapour pressure curves reported in the literature ex- 
tend to the range of desorption temperatures in these 
cxperimentsi2&2i2£ Due to the commonly observed in- 
crease of AH S with decreasing temperature we thus have 
to extrapolate and evaluate the vapour pressure curves 
at temperatures where desorption occurs in our experi- 
ments. This can be done by a fit to vapour pressure data 
using Antoine's expression*^ 
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FIG. 2: Linear fit to isothermal desorption data for naphtha- 
lene in the temperature range of 228 to 234 K. 

graphite. Table ^summarizes all parameters used for the 
computation of frequency factors as well as the tempera- 
ture at the desorption peak maximum used for computa- 
tion of the activation energy from the Redhead equation. 
The computed frequency factors are summarized in ta- 
ble [H] together with the resulting activation energies as 
obtained from eq. 



In(p(T)) = A 



D 



T + C 



(5) 



where A, B and C are the fit coefficients listed in tablcQ] 



TABLE I: Vapor pressures and desorption temperatures used 
to compute the frequency factor 



Molecule A B (K) C (K) 


o-(cm- 2 )° 


po(mbar) T max (K) 


Benzene 26 7640 30 


2.7xl0 18 


5.8x10" 


151 


Naphthalene 43 20f00 124 


1.6xl0 18 


1.3xl0 14 


235 


Coronene 37 30400 184 


0.9xl0 18 


2.6xl0 13 


390 


Ovalene 


0.6xl0 18 


7.6xl0 13 


490 



a See the text for references 

For ovalene, where no vapour pressure data is available, 
the pre-exponcntial factor 5.6 x 10 21 s _1 was obtained 
using the tabulated slope and offset from the Clausius 
Clapeyron equation in the form ln(p) = A — B/Tm^ 
The density a of benzene, naphthalene and coronene 
molecules adsorbed on the graphite surface can be ob- 
tained from LEED or STM dataiSiSl For ovalene where 
again no such data is available a can be estimated assum- 
ing that molecules are close packed with the polyaromatic 
rings oriented parallel to the surface. The resulting den- 
sities are summarized in table [I] 

The sticking coefficient s used in eq. iJSJ is here as- 
sumed to be close to unity which is commonly observed 
for homocpitaxial growth and should also be appropri- 
ate for adsorption of weakly interacting polyaromatics on 



2. Falconer- Madix analysis 

In this section we will perform an isothermal analysis of 
TD spectra of the type frequently referred to as Falconer- 
Madix method^* For this purpose we plot the logarithm 
of the desorption rate ln(—d6/dt) versus In 6* as evaluated 
at one specific temperature for several desorption traces 
of different initial coverage. A linear fit to the resulting 
data sets from sub-monolayer coverages as shown in Fig.[21 
can be analyzed according to: 



d6 

ln( ) =]nv + n\n6 

dt 



Eg 

k B T 



(6) 



If the intercept, / 



Eg 

k B T 



(7) 



is plotted versus 1 /T (see Fig.[3] for naphthalene) one 
obtains the pre-exponential frequency factor as well as 
the activation energy from a straight line fit to the 
data. From the slope of naphthalene and benzene des- 
orption isotherms we obtain the order of desorption n of 
0.95±0.02 and n=1.01±0.02 respectively. On the other 
hand, coronene and ovalene appear to follow fractional 
order kinetics with n=0.27±0.04 and 0.34±0.01 respec- 
tively. 

The prc-cxponcntials as well as the activation energies 
obtained from this analysis are summarized together with 
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FIG. 3: Plot of intercept I vs inverse temperature. 



the results from the Redhead peak-maximum analysis in 
Table CU 

Note, that the pre-exponential frequency factors in the 
present study are considerably larger than those com- 
monly used for desorption of smaller molecules where 
typical values are on the order of 10 13 s _1 to 10 15 s _1 . 
This reinforces the necessity to determine these parame- 
ters accurately and reliably either from experiments or by 



TABLE II: Binding energy and frequency factor : comparison 
between Redhead and Falconer-Madix analysis 



Redhead 
- 1 ) E a (meV) 



Falconer-Madix 
vis- 1 ) E a (meV) 



Benzene 



lxl0 16±3 0.50±0.08 5.0xl0 15±2 0.50±0.08 



Naphthalene 5x10 



0.8±0.1 l.OxlO 1 



0.90±0.07 



Coronene 
Ovalene 



2xl0 16±2 1.3±0.2 1.7xl0 18±0 - 5 1.5±0.1 



5xW 



2.2±0.2 



8.0xl0 18±0 - 5 1.97±0.08 



complementary theoretical investigations as recently re- 
ported by Fichthorn et al— The high pre-exponentials 
found here can qualitatively be accounted for by the 
larger difference between partition functions in the ad- 
sorbed and the transition state, if compared with des- 
orption of mono- or diatomic adsorbates, for example. 
A comparison of the pre-exponentials obtained here by 
the two different methods shows good agreement within 
estimated error bars. 



C. The interlayer cohesive energy of graphite 

The cohesive energy of a solid is commonly referred 
to as the energy needed to "disassemble it into its con- 
stituent parts'— while the work of cohesion is occasion- 
ally also referred to as the energy needed to "separate 
unit areas ..." of a medium "... from contact to infin- 
ity in vacuum" — The two are not identical and for a 
layered system with extremely anisotropic bonding like 
graphite the first can be identified with the exfoliation 
energy E ex , i.e. the energy needed to separate all layers 
of the crystal to infinity while the latter is equivalent to 
the cleavage energy E c i which is commonly slightly larger 
than the exfoliation energy. Previous work suggests that 
the energy needed to separate a single sheet of graphene 
from a graphite crystal, i.e. the exfoliation energy, is ap- 
proximately 18% smaller than the cleavage energy^ The 
normalization in this paper will be with respect to surface 
atoms and not to area. The area per surface atom in a 
graphene sheet is ((2.46-1CT 10 m) 2 - V3) = 1.05- lCT 19 m 2 . 

In the following we will use the results from the pre- 
vious section to determine the interlayer cohesive energy 
of graphite to which the dominant contribution is gen- 
erally believed to arise from long range van der Waals 
(vdW) interactions between graphene sheets. Our ap- 
proach for the determination of the interlayer cohesive 
energy - using the activation energies for desorption 
of small polyaromatic molecules, i.e. small 'flakes of 
graphene' - is based on the near-additivity of such vdW 
interactionsi22i2£ The latter is well established, for exam- 
ple, for the cohesive energy of alkanes where deviations 
from linearity in the number of chain segments are about 
1 % or less— The activation energy for desorption - as 
measured by TD spectroscopy - will here be identified 
with the binding energy of the adsorbate to the graphite 
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FIG. 4: Dependence of activation energy for desorption on the 
number of carbon atoms of four polyaromatic hydrocarbons. 



surface. The contribution of individual carbon atoms to 
this binding energy is derived from our data which - in 
the limit of infinitely large PAHs - would correspond to 
the energy needed to separate a graphene sheet from its 
parent crystal and is thus associated with the interlayer 
cohesive energy. Additional contributions to the desorp- 
tion energy from intcrmolccular interactions are assumed 
to be negligible. This is justified if adsorbate-adsorbate 
interactions are either small due to large mutual separa- 
tion - as in the case of the low coverage regime for ben- 
zene or naphthalene desorption - or this may be justified 
if adsorbate-adsorbate interactions within 2-D islands on 
the surface are comparatively small as in the case of ad- 
sorbed coronene or ovalene. 

A rough estimate of the exfoliation energy can be ob- 
tained by averaging over desorption energies per car- 
bon atom for all studied adsorbatcs which would yield 
67 meV/atom. In this case however, one neglects the 
small but significant contribution to the binding energy 
from hydrogen atoms and this approach is thus expected 
to overestimate the contribution of carbon atoms to the 
total binding energy. A more thorough analysis of the 
binding energies plotted in Fig.0] thus also has to ac- 
count for the contribution of hydrogen atoms. Here, this 
is done by optimizing the carbon-carbon and hydrogen- 
carbon interaction potentials to give optimum agreement 
of experimental energies with calculated binding ener- 
gies. Experimental data for this optimization are aver- 
aged over values obtained from both the peak maximum 
and isothermal analysis. The calculated binding energies 
are obtained in the usual manner by summation over em- 
pirical vdW pair-potentials: 



cc 



,|) (8) 



Lrn 



where the summation over n and / includes all carbon or 
hydrogen atoms in the adsorbate and the summation over 
m includes all atoms in the substrate. The carbon-carbon 
(C-C) and carbon-hydrogen (C-H) vdW potentials are 
those of the MM3 force field by Allinger and co- workers^ 
They have the form: 



U AB (r) = £ ab 



l 8 4000exp[-^V2.25f^) 6 
tab J \ r J 



(9) 

where for dissimilar atoms A and B cab is given by 
£ab = V £ A£B and tab is given by tab = ta + r B - 
The values (ec 5 ? 'c) an d (ch^h) given by Allinger and co- 
workers for carbon and hydrogen are (2.44 meV, 1.96 A) 
and (0.87 meV, 1.67 A), respectively. In the following 
we use the depth of the two vdW potentials V^ cc and 
V^ HC as free scalable parameters while the position of 
the potential minima is kept fixed. The binding energy 
Eg AK (Vf c , U HC , z) is then optimized by adjusting the 
molecule-surface distance z while the orientation of the 
aromatic rings of the molecules is fixed parallel to the 
graphite surface. We then determine the set of parame- 
ters V^ cc and V^ HC which gives best simultaneous agree- 
ment of all calculated with all experimental binding ener- 
gies by minimizing the corresponding root mean square 
deviation. Using the MM3 force field parameters the 
depth of the hydrogen vdW potential for interaction with 
a graphite surface is estimated to be 27meV/atom. If wc 
allow this to vary by at most ±5 meV we get best agree- 
ment with the experimental binding energies if the depth 
of the carbon-graphite potential is (52 ± 5)meV/atom. 
From Fig. 01 where calculated and experimental binding 
energies are compared one finds that deviations between 
experimental and calculated values are about 10% or less. 
The major uncertainty here arises due to the large error 
bars of the pre-exponential in the Redhead method and 
from the estimated error of the temperature calibration 
as well as the reproducibility of TD traces due to small 
variations of the sensed temperature for the Falconer- 
Madix analysis. For coronene and ovalene an additional 
error of about 2% arises due to uncertainties of the cov- 
erage calibration. The cleavage energy of graphite is ob- 
tained from the above value by accounting for the 18% 
higher energy previously observed for the separation of 
two crystal halves if compared with the separation of a 
single graphene layer from its parent crystal^ This yields 
a cleavage energy of (61 ± 5) meV/atom which is signifi- 
cantly larger than previously reported values. 

The earliest measurement of the exfoliation energy 
from heat of wetting experiments by Girifalco and Lad 
gave (260 ± 30) ergs/cm which - using the carbon atom 
density within a graphene sheet - is equivalent to (43 ± 
5) meV per surface atom£ Unfortunately, the only docu- 
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mentation for the heat of wetting data used for the deter- 
mination of E ex by Girifalco and Lad has been published 
in a thesis and is not readily available Also, no infor- 
mation on the kind of graphite sample used in the origi- 
nal experiments is available. Carbon powders commonly 
used for such experiments usually have higher surface ar- 
eas of the order of tens to hundreds of square meters 
per gram. They are usually treated with acids to oxi- 
dize non-graphitic contaminants and are then thermally 
treated to achieve the highest possible degree of graphiti- 
zation. The adsorptive properties and surface chemistry 
of such carbons, however, depend strongly on the history 
of the sample treatment^ An assessment of systematic 
and sample dependent uncertainties of the value reported 
by Girifalco is thus difficult. 

An even smaller value of (35l^g) meV/atom for the 
cohesive energy of graphite was obtained by Benedict 
et al. from the analysis of collapsed multi-wall carbon 
nanotubes^ This analysis is based on a measurement of 
the diameter of hollow 'bulbs' at the sides of three dif- 
ferent collapsed multi-wall carbon nanotubes with a pre- 
cision of about 1-2 A. N, other than statistical errors 
or those due to the limited accuracy of the bulb diame- 
ter determination may as well contribute to uncertainties 
associated with this value. This again makes it difficult 
to assess the relevance of possible systematic or statis- 
tical errors which could help to better understand the 
discrepancy between our and other experimental deter- 
minations of the cohesive energy of graphite. Note, that 
a distinct advantage of the present investigation is that 
experimental conditions and assumptions leading to the 
conclusions arc most clearly defined and that a well char- 
acterized model system is used for the investigation of the 
intcrlaycr cohesive energy of graphite. 

As stated in the introduction, theoretical estimates us- 
ing ab initio or semi-empirical methods show a much 
stronger variation from values as little as 8meV/atom to 
as much as 170meV/atomii*2iii2i£ili These large discrep- 
ancies are partly due to the inherent difficulties encoun- 
tered in the calculation of long-range dispersion forces. 
Even advanced calculations using non-local density func- 
tional theory which account for vdW interactions with 
reported values for a single pair of graphene sheets of only 
35 meV/atom& tend to underestimate the interlayer cohe- 
sive energies. However, the desorption energies reported 



here for different PAHs may serve as a useful benchmark 
for future studies to allow a better comparison of theo- 
retical binding energies with experiment. Reliable val- 
ues for the cleavage and interlayer cohesive energies may 
eventually be derived from a successful calculation of the 
interaction of PAHs with graphite if the interlayer forces 
in graphite can be treated on the same footing. 



IV. Summary 

In conclusion, we have presented a thermal desorp- 
tion study of benzene, naphthalene, coronene and ova- 
lene adsorbed the basal plane of graphite. Binding en- 
ergies were obtained by the peak-maximum method and 
alternatively by an isothermal analysis which also allowed 
to determine pre-exponential frequency factors. For the 
peak maximum method we derived pre-exponentials from 
vapor pressure data in combination with Antoincs law 
for extrapolation of the data to the temperature of des- 
orption. Both methods indicated that pre-exponential 
factors increase with adsorbate size from 10 14 s _1 to 
10 s _1 . The corresponding binding energies increase 
from 0.50 eV for benzene to 2.1 eV for ovalene. These val- 
ues were used to determine the interlayer contribution to 
the cohesive energy of graphite by assuming pairwise ad- 
ditivity of the interaction of carbon and hydrogen atoms 
to the total binding energy of the molecules. The re- 
sulting cleavage energy of graphite of (61 ± 5) meV/atom 
is derived from the average carbon atom contribution to 
the binding energy of the PAHs of (52 ± 5) meV/atom 
which can also be associated with the exfoliation energy. 
This is significantly larger than previous experimental de- 
terminations and these results provide an experimental 
benchmark for future theoretical investigations of inter- 
layer bonding and van dcr Waals interactions in graphitic 
systems. 
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